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Abstract 

O The flux of cosmic rays beyond the GZK cutoff (~ 10 20 eV) may be explained 

'-^ through their production by ultra high energy cosmic neutrinos, annihilating on the 



X 



relic neutrino background, in the vicinity of our galaxy. This process is mediated 



through the production of a Z boson at resonance, and is generally known as the 
Z-Burst mechanism. We show that a similar mechanism can also contribute to the 
super-GZK spectrum at even higher, ultra-GZK energies, where the particles produced 
at resonance are the Kaluza-Klein gravitons of weak scale mass and coupling from 
the Randall- Sundrum (RS) hierarchy model of localized gravity model. We call this 
mechanism Gravi-Burst. We discuss the parameter space of relevance to Gravi-Bursts, 
and comment on the possibility of its contribution to the present and future super- 
GZK cosmic ray data and place bounds on the RS model parameters. Under certain 
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assumptions about the energy spectrum of the primary neutrinos we find that cosmic 
ray data could be potentially as powerful as the LHC in probing the RS model. 
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1 Introduction 



About 25 years ago, Greisen, Zatsepin, and Kuzmin (GZK) noted that the observed spec- 
trum of proton, photon, and nucleus cosmic rays must virtually end at energies above ~ 10 20 
eV, the GZK cutoff [JEf. Their key observation was that Ultra High Energy Cosmic Rays 
(UHECR's) deplete their energy through various interactions with the 2.7° K Cosmic Mi- 
crowave Background Radiation (CMBR), over distances of order fO — f 00 Mpc. Above 10 19 
eV, nuclei are photo-dissociated by interactions with the CMBR, and a 10 20 eV proton loses 
most of its energy over a distance of ~ 50 Mpc. The analogous distance for a photon of the 
same energy is ~ 10 Mpc, due to e + e~ pair production on the radio background 

However, over the past three decades, different experiments have observed a total 
of about 20 events at or above this 10 20 eV bound[[|. Since there seem to be no feasible 
candidates for the sources of these cosmic rays, such as Active Galactic Nuclei, within a GZK 
distance ~ 50 Mpc of the earth, the observation of these events poses a dilemma. A number 
of proposals have been made to resolve this puzzle 0. One such proposal for the origin of 
the super-GZK events, due to Weiler, is based on the observation that UHECR neutrinos 
can travel over cosmological distances, with negligible energy loss||5], Therefore, if these 
neutrinos are present in the universe they could in principle produce Z bosons on resonance 
through annihilation on the relic neutrino background, within a GZK distance of the earth. 
The highly boosted subsequent decay products of the Z will then be observed as primaries 
at super-GZK energies, since they do not have to travel cosmological distances to reach us. 
This mechanism for producing super-GZK cosmic rays is referred to as Z-Burst. 

The Z-burst mechanism has the advantage that it does not assume physics beyond 
the Standard Model (SM) and is, therefore, minimalistic. However, any extension of the 
SM that provides a particle X which couples to vv and decays into the usual primaries can 
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in principle contribute to the super-GZK spectrum beyond the range presently observed. 
Assuming a mass m v ~ 10~ 2 — 10 _1 eV for neutrinos as suggested by atmospheric oscillation 
data, the particle X must have a mass of order the weak scale (~ 1 TeV) to be relevant to the 
spectrum near the GZK cutoff. In this paper, we will show that the massive Kaluza-Klein 
(KK) tower of gravitons in the Randall-Sundrum (RS) localized gravity model0] are viable 
candidates for particle X. 

The RS model is based on a truncated five-dimensional Anti-deSitter^ciSs) space- 
time, with two 4-d Minkowski boundaries. Our visible 4-d universe and all fields associated 
with the SM are assumed to be confined on one of these boundaries, referred to as the 
TeV brane, with the other 'Planck' brane boundary separated from us by a fixed distance 

r c ~ 10 M Pl , the compactification scale along the 5 th dimension; Mpi is the reduced Planck 
mass. The RS geometry is such that the induced metric on the visible TeV brane generates 
the weak scale from a 5-d scale M 5 ~Mpi, without fine-tuning, through an exponentiation. 
The interested reader is referred to Refs. |7], [8] for the details of the RS model and its numer- 
ous phenomenological implications. However, here we mention that a distinct feature of this 
model is that it predicts the existence of a tower of spin-2 KK gravitons, (n = 1, 2, 3, . . .), 
starting at the weak scale, and with weak scale mass splittings and couplings. 

Phenomenological studies || |9| suggest that the lowest lying KK graviton can 

be as light as ~ 400 GeV. The G^ have couplings to all particles, due to their gravitational 
origin and can be produced by vv annihilation, eventually decaying into qq, gg, 77, . . .. Thus, 
the G^ can in principle contribute to the super-GZK spectrum in a way that is similar to 
the Z-burst contribution. We call this graviton mediated process Gravi-Burst. 

Since the Z and G^ have different couplings and branching fractions to the observed 
primary particles, we expect that experiments may be able to distinguish between Z-burst 
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and gravi-burst initiated primaries. Also, depending on the behavior of the flux of neutrinos 
at super-GZK energies, more than one member of the KK graviton tower could contribute 
to gravi-burst. In this case, the RS model predicts a characteristic multi-peaked behavior 
for future data at super-GZK energies and beyond. However, collider experiments may be 
a better place to directly search for the graviton tower, with cosmic ray data providing 
complementary information as we will discuss in detail later. 

In the next section, we present the necessary formulae for estimating the super-GZK 
flux in the Z-burst model. We adapt this approach to gravi-burst and give the corresponding 
rate estimates in this scenario. Section 3 contains our results for a range of RS model 
parameters and a comparison with the Z-burst predictions. We will show that if the neutrino 
spectrum falls sufficiently slowly with energy we can use GZK data to greatly restrict the 
parameter space of the RS model. Our conclusions are given in section 4. 

2 The Burst Mechanism 

The burst mechanism relies on several well-motivated assumptions given the successes of 
the SM, Big Bang Nucleosynthesis and the observation of neutrino oscillations due to the 
existence of finite neutrino masses. This scenario is most easily demonstrated in terms of the 
conventional Z-burst. This model proposes that a high energy flux of neutrinos (and anti- 
neutrinos) are produced by some as yet unknown astrophysical source and collide with the 
relic background neutrinos in the galactic neighborhood. The origin of this flux is unspecified 
but constraints on its magnitude and energy dependence exist from Fly's Eye data|Hj. If 
the flux at the Z-pole is sufficient to explain the super-GZK excess then the Fly's Eye data 
tells us that the fall off with energy of the neutrino flux at somewhat lowers energies goes at 
least as fast at E~ 0,9 . A similar energy behavior may be expected above the Z-pole. 
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Due to the finiteness of neutrino masses one would expect that the local density of 
neutrinos will most likely be enhanced over the uniform cosmological background due to 
their gravitational clustering around the galaxyH. Massive neutrinos within a few Z widths 
of the right energy 

E? z = M 2 z /2m v = 4(0.1 eV/m„) x 10 22 eV , (1) 

will then resonantly annihilate into hadrons with the local anti-neutrinos (and vice versa) at 
the Z-pole with the large cross section 

< (J ann >^ = / TTn""'ami(s) = ^ 28nb , (2) 



M 2 anny ' 

where B% ~ 0.70 is the hadronic branching fraction of the Z. We assume that only left- 
handed neutrinos exist and employ the narrow width approximation. Given a neutrino mass 
hierarchy and the Super-Kamiokande atmospheric oscillation results [[FT]] we expect one of 
the neutrinos to have a mass near ~ 0.05 — 0.06 eV. (This follows from using the latest 
two parameter fit to the Super-K data which yields a value for Am 2 of 3.2 x 10~ 3 eV 2 
and by supposing that one of the neutrino masses is at least a few times larger than the 
second.) The locally produced 30 or so hadrons from the decay of the Z are then the effective 
primaries for the super-GZK events that are observed with energies in excess of ~ 10 20 eV. 
(In principle, there being three neutrinos, we should consider three different cases depending 
on their masses. This is a straightforward extension of the present discussion.) If the source 
of the initial neutrinos is randomly distributed in space then, as shown by Weiler||, we can 
calculate the total rate of super-GZK events induced by v — v annihilation at the Z pole 
within a distance D of the Earth as 

F z ~ E^F u {E Rz ) < a ann > z [° dx n{x) , (3) 

Jo 
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where the narrow width approximation has again been employed, F v (E Rz ) is the incident 
neutrino flux evaluated at the resonant energy, and n(x) is the column number density of neu- 
trinos. In deriving this expression it is assumed that the product < a ann > z J D dx n(x) << 1 
as is the case for the Z in the SM and in the RS model we consider below. In practice we are 
interested in rather close annihilation, i.e., values of D of order the GZK limit for protons 
which is ~ 50 Mpc. 

Weiler has shown that for reasonable ranges of the parameters the resulting value 
of the flux Fz can indeed explain the ~ 20 events beyond the GZK bound observed over 
the last few decades. We note that the model in its present form predicts that all of the 
super-GZK events are relatively well clustered in energy just beyond ~ 10 20 eV and that 
essentially no events should exist beyond those induced near the Z pole. Obviously, if such 
'ultra'-GZK events were observed then there must be new processes which can also lead to 
enhanced annihilation cross sections beyond those arising in the SM. 

In the RS model with the SM gauge and matter fields lying on the TeV brane there 
exist a Kaluza-Klein tower of massive, weak scale gravitons, G^ n \ with essentially electroweak 
couplings. There are basically two parameters in this model: the ratio c = k/Mpi, with k a 
mass parameter with a magnitude comparable to the five-dimensional Planck scale, and the 
mass of the lowest lying graviton state. The masses of the tower KK states relative to the 
first non-zero mode are given by the ratio of roots of the Bessel function J\ and c is expected 
to lie in the range ~ 0.01 — lj7|, ||, 9- Specifically, while the zero mode graviton couples with 

a strength M p \, all of the remaining KK tower states couple as A" 1 where A n = Mpi e~ 7TkTc . 
For values of kr c in the range 11-12, the RS model provides a solution to the hierarchy 
problem. The masses of the KK states, G^ n \ are then given by m n = kx n e~ nkVc with x n 
being the n th roots of J\. This implies that the tower mass spectrum is completely determined 
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once the mass of the lowest lying excitation is known and is given by m n = m\x n jx\. Thus 
we see that the parameters c = k/ M Pi and mi determine all of the other quantities within 
the RS model. 

Both phenomenological and theoretical constraints can be used to restrict this two 
dimensional model parameter space as has been discussed in our previous works [^, [jj. As 
an initial numerical example of the Gravi-burst mechanism let us consider the specific case 
where mi = 600 GeV and k/Mpi = 0.1 which is a point in parameter space that is allowed 
by all of the current constraints; it is then straightforward to consider the more general 
case. For these specific values of the parameters we can calculate the cross section for vv 
annihilation into hadrons. This cross section now has a number of distinct contributions 
besides those arising from qq final states as in the decay of the Z in the SM. (Even in the 
SM, above threshold, neutrinos can annihilate into pairs of W's, Z's and top quarks which 
can subsequently hydronically decay.) In the RS case, the gravitons not only lead to qq, tt, 
W + W~ and ZZ final states but also to pairs of gluons and Higgs bosons, gg, hh. Gluons 
fragment directly into hadrons while the SM Higgs bosons decay mostly to bb. (Here we 
assume for numerical purposes that the mass of the Higgs is 120 GeV.) By combining all of 
these individual process cross sections, including interference with SM Z exchanges, we can 
calculate the full energy dependence of the total vv — > hadrons cross section. This allows us 
to determine the ratio of expected cosmic ray rates for super- and ultra-GZK events in units 
of the Z-pole induced rate Fz computed above. It is important to note that in forming this 
ratio almost all of the astrophysical uncertainties cancel except for the energy dependence 
of the neutrino flux. We find 

D , ^ _ Fsm + grav(V~s) _ 2V~s cr! Z +GRAV (VS) (Mz/VSy 

Ml < a ann >z ' 1 ) 

where we have assumed that the neutrino spectrum above the the resonant Z pole energy 
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falls in a power-like manner as ~ E~ x / 2 . We denote by Fsm the complete energy dependent 
flux anticipated in the SM beyond that obtained through the use of the narrow-width ap- 
proximation alone. (In what follows it will be sufficient to assume that this power-like fall off 
adequately describes the neutrino spectrum for a few orders of magnitude in energy above 
E^ z .) Integration of R over a range of ^/s values then tells us the relative rate of events 
expected in the RS model to those originating from Z-bursts. 
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Figure 1: Energy weighted total cross section for hadron production in units of that for the 
Z pole in the Weiler Z-burst model for A = as a function of center of mass energy for 
the SM(the relatively flat lower blue curve) and in the RS model(the upper red curve) with 
c = k/Mpi = 0.1 and m\ = 600 GeV. The small irregularities in the curves are due to WW, 
ZZ, hh and tt thresholds. 



To get an idea of what this ratio looks like as a function of energy we show the simplest 
specific case where A = in Fig. 1. Note that the integral of R under the Z pole gives the 
value unity as it should to reproduce the Weiler results. We are also reminded by the figure 
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that even in the SM there exists a long high-energy tail to this ratio. In the RS scenario, the 
Z peak is followed by a number of ever widening graviton peaks which also yield reasonably 
large cross sections. From the figure, one can see that if the neutrino flux falls off slowly 
enough with energy we should expect events at even higher energies than those observed at 
present assuming that the Z accounts for the 'usual' super-GZK events. We will make this 
assumption in what follows, i.e., that the Z-burst scenario explains the observed super-GZK 
events. Given that hadronic multiplicities grow only very slowly with y/s, as is observed 
in e + e~ annihilation data, we would interpret events induced by Gravi-bursting on the first 
graviton resonance to result in hadronic effective primaries that have energies approaching 
10 22 eV. As of yet no such events have been recorded which places bounds on the allowed 
parameters of the RS model for different values of the neutrino energy spectrum described 
by the parameter A. 



3 Analysis 

As we found in the last section, under the assumption that Z-bursts explain the super-GZK 
events the existence of even higher energy ultra-GZK events is a rather generic prediction 
of the RS model. Let us restrict ourselves to the region y/s > 300 GeV which corresponds 
to effective primary energies in excess of 10 21 eV of which none have been yet observed. 
Integrating R above this lower bound, even in the SM, can yield some 'background' events; 
the use of the narrow width approximation is not strictly correct in that some rare events 
can arise from values of y/s away from the Z pole. In the RS case, we integrate R over the 
region from 300 GeV up to y/s = 4mi beyond which perturbation theory fails. f\ This yields 

a conservative lower bound on the total number of ultra-GZK events that are predicted in 

^We note that the RS model as described in four dimensions is a non-renormalizable theory. In addition, 
once the value of y/s significantly exceeds A w the theory also becomes non-perturbative and only qualitative 
statements can be made about the behavior of the cross section g. 
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the RS model since most certainly more events can arise from even larger values of ^/s. 
Integrating R in the SM over the above ranges and assuming that the 20 super-GZK events 
are from the Z-pole region, we find that for A = 1(2, 3) we would expect to have already seen 
~ 0.24(0.04, 0.008) ultra-GZK background events from the tail of the Z pole, which is quite 
acceptable. (As discussed above we might expect that A > 1.8 is allowed by Fly's Eye data 
if the energy dependence of the neutrino spectrum below and above the Z-pole are similar.) 
Performing the same calculation in the RS model for a fixed set of values of m\ and A it will 
be clear that for some range of k/Mpi the cross sections will be too large to have avoided 
the present non-observation of ultra-GZK events. In the usual manner this means that we 
can place a 95% CL bound on k/Mpi as a function of mi for different assumed values of A, 
using the existing data. 

The dashed curves in Fig. 2 show the results of this analysis using the existing data 
for various values of A; other constraints obtained on the (k/Mpi)-mi parameter space from 
our earlier work|J are also shown. We see immediately that the effectiveness of the bound 
is quite strongly dependent on the value of A. For A > 3 at most only a tiny region beyond 
that excluded by existing Tevatron Run I data is now ruled out. As A decreases the size 
of the region presently excluded by cosmic ray data grows rapidly. For A = 2 a substantial 
region allowed by the present Tevatron data becomes excluded. Furthermore, a sizeable 
region beyond that accessible at Run II with a luminosity of 30 fb~ l is also excluded. Using 
accelerators alone this region would be inaccessible until after the LHC turns on but here 
we see that it would be excluded by cosmic ray data provided A < 2. For A < 1 the bound 
is extremely powerful and at most only a tiny sliver of the RS model parameter space would 
remain viable. 

In the region below the dashed curves, which is not excluded by existing cosmic 
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Figure 2: Allowed region in the {k/M P j)-mi plane. The solid(dotted) diagonal red curve 
excludes the region above and to the left from direct searches for graviton resonances at 
the Run 1(11, 30 fb^ 1 ) Tevatron. The light blue(green) curve is an indirect bound from 
the oblique parameter analysis (based on the hierarchy requirement that < 10 TeV) 
and excludes the region below it. The black dashed(dotted) curves excluding the regions 
above them at 95% CL based on present (anticipated future Auger) cosmic ray data. The 
top(bottom left, bottom right) panel corresponds to A = 3(2, 1) which describes the fall with 
energy of the neutrino flux as E~ x l 2 . 
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ray data, we might expect ultra-GZK events to show up in future experiments at reasonable 
rates. If this does not happen it's clear that the present bounds discussed above will improve 
drastically especially with the new cosmic ray observatories such as Auger [O coming on 



line. Within a 5 year period of data taking at Auger one would expect ~ 1000 super-GZK 



event s ||13[] induced by Z-Bursts with correspondingly higher sensitivity to the ultra-GZK 
region. If no events above the SM background from the Z pole tail are observed at Auger 
during this period we can repeat the analysis above to obtain strong constraints on the 
RS parameter space as shown by the dotted curves in Fig. 2. We find the SM background 
expectations in this case for A = 1(2, 3) to be ~ 12.5(2.05, 0.422) events. Here we see that for 
A = 2, 3 the size of the presently allowed region is quite significantly reduced. Particularly 
note the case A < 1 where we find that the non-observation of any ultra-GZK events at 
Auger would completely exclude the RS model with the SM gauge and matter fields on the 
wall. This is a very powerful result. 

If events above background are observed at Auger due to gravi-bursts they will have 
two distinctive characteristics. First, due to the resonance structure predicted by the RS 
model the energies of the effective hadronic primaries will show peaking at a set of fixed 
energies provided the energy resolution of the detectors is sufficiently good. Second, in 
addition to the rather 'soft' photons arising from conventional fragmentation 7r°'s, much 
harder photons can arise from the direct decays of the gravitons in the KK tower. As shown 
in our earlier work, gravitons in the mass range of interest can decay with a reasonable 
branching fraction, ~ 4 — 5%, into photons. Since they carry half of the energy of the 
resonance mass these photons will have energies an order of magnitude or more larger than 
those arising from 7r° decays. Of course with this rather small branching fraction a reasonable 
number of more 'ordinary' ultra-GZK events should be observed before one induced by these 
very hard photons. 



12 



So far we have only discussed the case of a single massive neutrino; data based on 



oscillation solutions to the solar neutrino problem [fUJI suggest a second massive state exists 
near 3 x 10~ 3 eV, about a factor 20 or so in mass below the 0.06 eV case discussed above. 
This second neutrino can also induce a Z-burst but only if the energy of the corresponding 
incident neutrino is ~ 20 times larger, ~ 1.4 x 10 24 eV. In comparison to that for the case 
of the higher mass neutrino the flux of these lighter neutrinos would be ~ (20) A / 2 times 
smaller. This assumed neutrino mass ratio, strongly suggests that A > 2 to avoid ultra-GZK 
events from the second Z resonance. To eliminate any additional background generated by 
this new Z contribution, we would need to raise lower bound on our y/s integration, which 
is actually an integral over the neutrino energy. (The \/s = 300 GeV lower bound translates 
into a minimum neutrino energy of ~ 8 x 10 23 eV and thus would now include additional 
background events from the second Z peak.) Raising the minimum neutrino energy by a 
factor of two would remove the second Z contribution while still staying comfortably below 
the excitation energies of any of the gravitons. In this case we would expect at most only 
slight alterations in the bounds presented above. 

Before we conclude we briefly discuss a generalization of the RS model where the SM 
gauge and matter fields are taken off the TeV brane|| and how this would influence our 
results above. In this case not only can graviton towers be exchanged in the vv — > hadrons 
process but now there are also Z boson towers whose members are generally interspaced in 
mass with the gravitons. If these additional contributions are also present it is quite possible 
that the neutrino cross section can be significantly enhanced leading to even stronger limits 
than those obtained above using current data. One might also expected that with increased 
cross sections it might be possible to probe cases where the slope of the neutrino energy 
spectrum is even steeper than what we have considered here. Unfortunately, to determine 
how much our previous results are modified in a quantitative manner requires a detailed 
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analysis which is far beyond the scope of this paper. 

4 Conclusions 

In this paper we have examined the possible contribution to the spectrum of cosmic rays 
beyond the GZK cutoff due to new physics arising in the Randall-Sundrum model of localized 
gravity. Our analysis is based on the assumptions (i) that the events observed immediately 
above the GZK bound can be explained by the Z-Burst mechanism and (ii) the neutrino 
spectrum needed for Z-Bursts extends s few orders of magnitude further in neutrino energy 
with a reasonably slow fall-off. If these conditions hold then the existence of a series of s- 
channel Kaluza-Klein graviton resonances in the vv — > hadrons channel, which is predicted in 
the RS model, can lead to events with even higher energies, ultra-GZK, due to Gravi-Bursts. 
The rate for these bursts are generally at or near the present level of observability for a wide 
range of RS model parameters. The fact that such events are not as yet observed can be used 
to constrain the parameter space of the RS model once a specific form of the neutrino energy 
spectrum is assumed. These bounds can be more restrictive than those that can be obtained 
from the lack of graviton resonance production at the Tevatron during Run II (30 fb" 1 ) if 
the fall-off with energy of the UHECR neutrino flux is linear or less steep. If ultra-GZK 
events are not observed by future experiments such as the Auger Array, then the resulting 
bounds on the RS model can be complementary to those obtainable at the LHC. If such 
events are observed at future experiments, the RS resonance structure may be observable 
given both sufficient statistics and good hadronic energy resolution. In addition to hadronic 
modes, the RS graviton KK tower states can directly decay to photon pairs which will have 
more than an order of magnitude greater energies that those that can arise due to ordinary 
fragmentation into 7r°'s which subsequently decay into two photons. If photon and hadron 
induced showers can be distinguished at such energies this will provide a unique signature 
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for the RS model as the origin of the ultra-GZK events. 
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